1 



Phase Transitions in Lattice QED 

M. Baig" 

"Grup de Fi'sica Teorica, IFAE 
Universitat Autonoma de Barcelona 
08193 Bellaterra (Barcelona) SPAIN 

The main objective of the work presented here is to understand the appearance of phase transitions in pure 
gauge and scalar lattice QED, both in the compact as well as the noncompact formulations. 



1. INTRODUCTION 

Recent developments in the study of lattice 
QED without fermions have lead to a rather puz- 
zling scenario 

• Pure gauge compact QED with periodic 
boundary conditions (PBC) shows a first 
order phase transition associated to a 
monopole condensation phenomena [1,2]. 

• Pure gauge non compact QED with PBC 
is a gaussian model that does not exhibits 
phase transitions, but shows monopole per- 
colation [3]. 

• Non-compact QED with PBC coupled to 
unitary Higgs fields shows two different 
lines, one of second order phase transitions 
and another related to a monopole percola- 
tion phenomena [4]. 

We will present here four topics that are in- 
tended to clarify some of the relations between 
monopoles and phase transitions. 

2. MONOPOLE PERCOLATION 

The lattice action used in our simulation [5] is 
the original Wilson (compact) action 

Sgauge ='^COsQ^^{n). (1) 

n ^ly 

We have applied the standard Metropolis algo- 
rithm to simulate the compact U(l) lattice gauge 
theory with standard periodic boundary condi- 
tions on lattices from S'* to 14'*. We have made 



measurements of the monopole observables near 
the phase transition. 

The evolution of the internal energy, the cluster 
susceptibility x and the rimax/ntot order parame- 
ter in "Monte Carlo time" , i.e. the iteration num- 
ber, show a clear correlation. The "tunneling" 
between the two phases is very clear for all three 
observables at the same point. This phenomena 
is characteristic of a first order phase transition. 

Our conclusion is that in compact pure gauge 
QED, the percolation threshold occurs at the 
same point as the deconfining phase transition, 
i.e. the point of monopole condensation. More- 
over, the Umax I ntot order parameter used to char- 
acterize monopole percolation proves to be an ex- 
cellent order parameter to study the location of 
the apparent first order phase transition. 

3. FIXED BOUNDARY CONDITIONS 

Numerical simulations of lattice gauge theo- 
ries are actually performed in a finite space-time 
box, and computer limitations imply severe re- 
strictions on the lattice size of this box. To avoid 
the border "effects" generated by such a box, pe- 
riodic boundary conditions are usually adopted. 
From a geometrical point of view, this fact im- 
plies that the system is actually simulated on an 
hypertorus. 

A possible explanation for the observed discon- 
tinuity in the internal energy on pure gauge com- 
pact QED with PBC relies on the fact that topo- 
logically nontrivial loops which wrap around the 
lattice are permitted. Remember that in _D = 3 
monopoles are 0-dimensional pointlike excitations 



2 



whilst in D = 4 they become 1-dimensional, and, 
due to the magnetic flux conservation, they form 
closed loops. 

In order to avoid toroidal topology we have im- 
posed to our lattice fixed boundary conditions 
(FBC). This means that we fix -and keep- to the 
unity all links belonging to the boundary of the 
lattice. From a geometrical point of view this is 
equivalent to a collapse to a single point all the 
boundary obtaining an effective spherical topol- 
ogy. 

Our numerical results [6] come from some 
thermal cycles performed for lattice sizes 
6^, 8^, 10^, 12^, and 16^. No sign of discon- 
tinuity has been observed. In addition, the emer- 
gence of a peak of the specific heat which grows 
with L and moves towards /3 = 1 has been clearly 
observed. Also the growing of that peak seems 
slow enough to ensure the absence of a weak first 
order transition. Furthermore, the Binder cumu- 
lant shows a clear convergence to the expected 
value of 2/3 characteristic of a second order tran- 
sition. 
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Figure 1. Binder cumulant measurement. 



These results, although being basically quali- 
tative, point out the continuous character of the 
U(l) phase transition when periodic boundary 
conditions and the spurious effects they bear with 
are avoided. 



4. NON-COMPACT SCALAR QED 

The non-compact coupled gauge-Higgs model 
(scalar QED) has been extensively studied in a 
recent paper [4]. Our model has gauge fields in 
the links of the lattice that are elements of the 
U{1) gauge group O^^fi and matter fields on the 
sites (f)x = e'"'-^-'. The action is 

^ = ^/5E^p-^E('^^^-.M'^-+M + cc), (2) 

where 6p is the circulation of the gauge field 
around a plaquette, /3 = i is the gauge cou- 
pling and A the Higgs coupling. Phase diagram 
shows two separate regions, a Coulomb phase and 
a Higgs phase separated by a second order transi- 
tion line. In the limit 7 — cxd limit the transition, 
however, becomes of first order. This limit can be 
interpreted as an Integer Gaussian Model (IGM) 
since the action becomes 

^oo = -f E(2™p)'- (3) 

p 

It is interesting to note that eq. 3 is just the 
Villain action for the pure gauge compact QED 
but with a 27r^/3 instead of l/l3y. Moreover, if 
one consider the Lagrangian version of the Hamil- 
tonian Loop Representation for the U(l) theory, 
one obtains an IGM dual eq. 3. This fact im- 
plies a new view of the nature of the phase tran- 
sition of pure gauge compact QED. A numerical 
simulation has been performed [7] observing an 
apparent first order phase transition, as expected 
since this model has the same topological struc- 
ture than pure gauge compact QED. 

5. FIELD THEORIES OF THE THIRD- 
RANK 

Some years ago, Poliakov pointed out that an 
abelian gauge theory with fields of the third rank 
will present, in four dimensions, point-like topo- 
logical excitations and hence, no phase transitions 
coming from monopole condensation. 

As a first step in the analysis of this question 
we have studied the four-dimensional abelian Z(2) 
case [8]. We denote by n a lattice site and by a 
link starting from point n in the direction /j, of an 
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4-d hypercubical lattice. Starting from this point 
there are 4 positive links. A simple plaquette de- 
fined by directions /j,, i> is denoted by n^j^. We 
assign a gauge variable (ctj- = ±1) to each pla- 
quette. Elementary three-dimensional cubes are 
composed by six plaquettes (i G dcube, i-e. the 
plaquette belongs to the perimeter of the cube). 
Three of them are the plaquettes associated to the 
combinations of the three index /U, J^, ?y related to 
the links starting from a given point n in the pos- 
itive directions. The remainder three plaquettes 
are those that "close" the cube. The partition 
function is defined, then, as: 

Z=Y,e-P^- E= J2 E{cuhe)- (4) 

{ctj} cu6e = l 

where 

E(cuhe) = \- W cji- <Ti=±l, (5) 

with /3 > being the inverse temperature in nat- 
ural units. 

A duality transformation [9] relates this four 
dimensional gauge model to the four dimensional 
Ising model, which has a second order phase tran- 
sition with mean field critical exponents [10]. 

Our estimates for the critical exponents v and 
a exclude a first order phase transition, and they 
are consistent with those expected from the 4-d 
Ising model. However, a plot of the energy his- 
togramjust over the phase transition point, shows 
a clear two-peaks structure. Moreover, a further 
finite-size analysis of the p™^\ reveals that in 
the limit L ^ co the surface tension is compatible 
with zero, giving place to a continuous transition. 

We conclude that the system actually experi- 
ments a second order phase transition with mean 
field critical exponents, as predicted by duality, 
and that the intriguing double peak structure 
seems to be due to an added phenomena, that 
may be understood in terms of the topological ex- 
citations of the model, i.e. the monopoles. This 
image is consistent with the fact than in five di- 
mensions, where the topological excitations be- 
come 1-dimensional, this model exhibits a clear 
first order phase transition. 



6. CONCLUSIONS 

Our main conclusions can be summarized as 
follows 

• Pure gauge compact QED with PBC shows 
a monopole percolation phenomena coupled 
to the monopole condensation. 

• Pure gauge compact QED with Fixed 
Boundary Conditions (spherical topology) 
shows a behaviour compatible with a sec- 
ond order phase transition 

• The 7 — cxD limit of Non Compact Scalar 
QED with PBC is related to the Villain 
form of compact QED and also to the Loop 
Model for QED, showing a first order phase 
transition. 

• Numerical a analysis of a third rank abelian 
Z(2) theory shows an intriguing behaviour 
over the phase transition that may be orig- 
inated by point-like topological excitations. 

The simulations done here used mainly the 
CRAY C90 at PSC and the CRAY Y/MP of 
CESCA. 
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